The WiggleZ Dark Energy Survey: constraining the evolution of Newton's constant 

using the growth rate of structure 
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We constrain the evolution of Newton's constant using the growth rate of large-scale structure 
measured by the WiggleZ Dark Energy Survey in the redshift range 0.1 < z < 0.9. We use this 
data in two ways. Firstly we constrain the matter density of the Universe, fim (assuming General 
Relativity), and use this to construct a diagnostic to detect the presence of an evolving Newton's 
constant. Secondly we directly measure the evolution of Newton's constant, GcH, that appears in 
Modified Gravity theories, without assuming General Relativity to be true. The novelty of these 
approaches are that, contrary to other methods, they do not require knowledge of the expansion 
history of the Universe, II{z), making them model independent tests. Our constraints for the second 
derivative of Newton's constant at the present day, assuming it is slowly evolving as suggested by Big 
Bang Nucleosynthesis constraints, using the WiggleZ data is Gcff (io) ~ —1.19 ± 0.95 ■ 10~^''ft^yr~^, 
where h is defined via Ho = 100 h km s~^ Mpc~^ , while using both the WiggleZ and the Sloan 
Digital Sky Survey Luminous Red Galaxy (SDSS LRG) data is Gcff(io) = -3.6 ±6.8 • 10"^^/i^yr"^, 
both being consistent with General Relativity. Finally, our constraint for the rms mass fluctuation 
as using the WiggleZ data is as = 0.75 ± 0.08, while using both the WiggleZ and the SDSS LRG 
data as = 0.77 ± 0.07, both in good agreement with the latest measurements from the Cosmic 
Microwave Background radiation. 



I. INTRODUCTION 



Recent observations suggest that the Universe is un- 
dergoing a phase of accelerated expansion, usually at- 
tributed to an unknown ideal fluid dubbed dark energy. 
Most of the evidence for the existence of dark energy 
comes from geometric tests that measure the expansion 
rate of the universe H(a) = a/a at different epochs, 
where a{t) is the scale factor of a Friedmann-Robertson- 
Walker metric. Examples of such tests include measure- 
ments of the luminosity distance, di(a), using standard 
candles like Type la supernovae (SN la) ^1), ^ and mea- 
surements of the angular diameter distance, (Ia (a) , using 
standard rulers such as the scale of the sound horizon 
at last scattering [3] and baryon acoustic oscillations 4] . 
Even though these tests are presently the most accurate 
probes of dark energy, the mere determination of the ex- 
pansion rate, H{a), is not able to provide significant in- 
sight into the properties of dark energy to distinguish it 
from models that attribute the accelerating expansion to 
modifications of general relativity [51-Q . 

A smoking gun signature of modified gravity theories, 
such as f{R) gravity [9] or scalar-tensor models 10], is 
the fact that they predict that Newton's constant Gm 
is evolving with time. Generically, modified gravity the- 
ories can be cast in such a way as to have an effective 
Newton's constant that changes with time Goff(i) [lll - 
|l3|. However, this variation cannot be arbitrarily large 
and has already been constrained by Big Bang Nucle- 



osynthesis to be within 10%. See, for example, Ref. |14| 

which gives G'cff(tsBA')/Goff (to) = 1.09±[5:?i, where ^bbn 
is the time of Big Bang Nucleosynthesis, while io is the 
present day. This allows us to make the simplifying as- 
sumption that Geff is a slowly varying function. Fur- 
thermore, Solar System tests and SN la observations 



(see [15| and references therein) suggest that the first 
time derivative Gog (to) is almost zero and that for the 
second time derivative we have \g2\ ^ G(l) , where 

52=Gcff(io)/(Geff(io)i?o')- 

The additional observational input that is required to 
discriminate between dark energy and modified gravity 
theories is the growth function of the linear matter den- 
sity contrast (5 = — , where p represents the background 
matter density and 5p its first order perturbation. The 
reason for this is the fact that the growth of cosmic struc- 
ture is a result of the motion of matter and therefore is 
sensitive to both the expansion of the Universe H{a) and 
Gcff. 

More specifically, it can be shown that in modified 
gravity theories on sub-horizon scales, ie k/^ ^ a^H^ 
where k is the wave-number of the modes of the per- 
turbations in Fourier space, the grow th factor satisfies 
the following differential equation [lll - llSl [l6| : 
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where primes denote differentiation with respect to the 
scale factor, H{a) = - is the Hubble parameter and we 



assume the initial conditions i5(0) = and 6'{0) = 1 for 
the growing mode. When Geff(a) = 1 we get GR as 
a subcase, while in general in modified gravity theories 
Goff(a) is time and scale dependent. For example, in 
f{R) theories we have that \l^ 
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which reduces to GR only when f{R) = R — 2A. 

The exact solution of Eq. (II. ip for a flat GR model 
with a constant dark energy equation of state w is given 
for the growing mode by [13, [l3 : 
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for H(af/H^=nma-^ + {l-n,^)a-^('^+"'\ (1.5) 



where 2Fi{a, b; c; z) is a hypergeometric function defined 
by the series 
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on the disk \z\ < 1 and by analytic continuation else- 
where, see Ref. [1^ for more details. 

In more general cases it is not possible to find an ana- 
lytical solution for Eq. (jl.ip . but in Ref. 20] it was shown 
that the growth rate /(a) = j^ can be approximated 
as 
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where this approximation for 7 is valid at first order for 
a dark energy model with a constant equation of state w 
and for ACDM we have 7 = ^^ ~ 0.55. Conveniently, the 
approximation of Eq. (jl.71) can be used to accurately fit 
a wide variety of different scenarios, including modified 
gravity models, if 7 is allowed to be a free parameter. 

In order to get a prediction for GcS from the growth 
rate data one would need to solve Eq. (|l.ip and this can 
only be done for a specific theory, as knowledge of -ff (a) 
is required. Therefore, in order to reconstruct GcS the 
usual procedure would be to assume a model for H(a), 
fit the approximation of Eq. (|1.7p to the data and finally 
use Eq. (|l.ip to get a prediction for Gcff . This obviously 
has several drawbacks, as the growth data alone are not 
sufficient to get viable fits to H{a) and are usually com- 
plemented by SN la and CMB data. The SN la data 



actually determine the luminosity distance ^^(a), there- 
fore converting these measurements to H(a) depends on 
differentiation of noisy data and an assumption about the 
spatial curvature VIk and this obviously results in large 
error bars for the best fit parameters. Finally, any result 
would ultimately depend on the specific choice of H{a), 
or in other words it would suffer from model-choice bias, 
whereas what we propose is a model independent tech- 
nique that can be used even without the CMB and SN la 
data. In order to solve these problems we constructed a 
diagnostic that is based solely on data, with no assump- 
tions about the expansion history of the Universe, that 
is able to detect the presence of an evolving Newton's 
constant. We describe both our new method and the 
diagnostic in detail in Section Hill 



II. THE WIGGLEZ DATA 

The WiggleZ Dark Energy Survey at the Anglo- 
Australian Telescope [21| is a large-scale galaxy redshift 
survey of bright emission-line galaxies mapping a cos- 
mic volume of order 1 Gpc^ over the redshift interval 
z < 1. The survey has obtained of order 200,000 rcdshifts 
for UV-selected galaxies covering of order 1000 deg^ of 
equatorial sky. In this paper we analyze the subset of 
the WiggleZ sample assembled up to the end of the lOA 
semester (May 2010). We include data from six survey 
regions in the redshift range 0.1 < z < 0.9 - the 9-hr, 11- 
hr, 15-hr, 22-hr, 1-hr and 3-hr regions - which together 
constitute a total sample of A^ = 152,117 galaxies. 

In Ref. 22] the authors presented fits for the growth 
rate of structure within this sample using redshift-space 
distortions in the 2D power spectrum. These fits in- 
cluded a full exploration of the systematic errors arising 
from the assumption of redshift-space distortion models 
based on perturbation theory techniques, fitting formulae 
calibrated by N-body simulations, and empirical mod- 
els. The WiggleZ data were analyzed in four redshift 
slices with effective redshifts z = (0.22,0.41,0.6,0.78). 
The growth rate measurements determine the quantity 
f{z) as{z) = (0.42 ± 0.07, 0.45 ± 0.04, 0.43 ± 0.04, 0.38 ± 
0.04) at these four redshifts. The reader is referred to 
Ref. [22 for full details of how these measurements were 
performed. 



III. THE METHOD 

Instead of solving Eq. (|l.ip for S{a), we will solve it for 
H{a) [Hill]. Then we get ^ 
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^ See Appendix |X] for more details. 



The growth rate is defined as /(a) = ^j^ and therefore 
S can be written in terms of / as (5 = (5(l)eJi /(^)/^<*^. 
However, WiggleZ gives the data in the form f -as, where 
the rms mass fluctuation as{z) can be expressed as [20| 
as(z) — '?A s S{z) and as a result we can express the 

product as /(a) • as{a) — -§(jy a S'{a). So, it would be 
convenient to define a new parameter 

fs8ia) = f{a)-as{a)^^aS'{a) (3.2) 

where f = ^frjy- Using this we get for the growth factor 
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while using Eqs. (|3.ip and (|3.3p we get for the Hubble 
parameter H{a) 

a*fs8{aY Jq Jo y 

(3.4) 
The parameter <^ has canceled out, so evaluating Eq. p.4p 
at a = 1 and solving for D,^ we get 
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where 



/s8(l) Jo y/s8(l) 

Notice that Eq. p.Sp does not include the Hubble param- 
eter H{z) or any assumption of some dark energy model 
at all. Therefore, the calculation of ilm or Gcfi can be 
carried out by using only the growth rate data. Now, 
Taylor expanding GcS around a — 1 
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where we have set f/„ = — da"^ |a=i; and plugging it 
into Eq. p.Sp we get 
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where 
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Clearly, the constants /„ depend solely on the data. We 
can get another constraint by demanding that Gcff(a = 
0) = 1, then using the Taylor expansion we get 
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(3.10) 



If we don't make the assumption that Geff (a ^1)^1 
then it can be shown that the growth of perturbations 
(5(a) in a matter dominated Universe (J7m = 1), but with 
a non-constant Geff, does not scale like S{a) ^ a for a <C 1 
as expected^. Also, since modified gravity models are 
used as alternative explanations of dark energy, which is 
prevalent at late times, any deviations from GR would 
also be expected at late times not at early times. We can 
use Eq. p.lOp to eliminate go, since 
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and from Eq. (|3.8p we get 
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or by combining these two and rearranging 
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In GR we have go = \ and (?„ — for n > 1, so we can 
define the GR value of fim as 
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(3.14) 



while if we could determine the form of Goff (a) accurately 
then we would define the "real" value of ilm as 
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So, if we form the ratio between the two parameters we 
have 
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(3.16) 



^ See Appendix [B] for analytic formulas for S{a) in a matter dom- 
inated Universe with a non-constant G^s- 



Also, Eq. (|3.14p can be used to construct a diagnostic 
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(3.17) 



where fim is the vahie of the matter density as measured 
from other independent observations. Therefore, if we 
assume that the value of Om is independently and accu- 
rately determined by other observations, then any statis- 
tically significant deviation of the quantity from zero 
or the ratio r from 1, clearly and uniquely identifies an 
evolving Gofr and consequently favours modified gravity 
theories. Finally, notice that the diagnostic depends 
solely on measurable quantities and not on any model. 

For concreteness we will also consider a specific model 
for Gcff. However, any variation of Newton's constant 
can be constrained by BBN to be less than 10%, see 
for example Rcf. JJJ who report Gcff(tBBN)/Geft(io) — 
1.09±Q-^g, which motivates us to assume that Gcs is a 
slowly varying function. Also, the Solar System tests 
and SN la observations (see jl5| and references therein) 
suggest that I51I ~ lO'^h'^ < 1 and \g2\ ^ 0(1), 
so we can safely neglect gi. However, the measurable 
quantities are Gcff (io) and Gcs{to), so we should express 
these (after a bit of algebra) as Goff(to) = 9iHq and 
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52 ) • So, after neglecting 
gi , the expression for the second derivative can be simpli- 
fied to Gcff (to) — 52-ffo- Furthermore, modified gravity 
models are used as alternative explanations of dark en- 
ergy, which is prevalent at late times, so any deviations 
from GR would be expected at late times not early times. 
As a result we will demand that Gcff (a ^ 1) ^- 1. 

Additionally, we will not consider models with wildly 
oscillating Gcff. The reason for this is that since Gcff is 
calculated via differentiations of the gravity Lagrangian, 
for example in the /(i?) case see Eqs. (|1.2|) - (|1.4p . 
then one would need to have a Lagrangian with terms 
that exhibit explicit oscillatory behavior. In the case 
of the f{R) theories one could consider the Lagrangian 
f{R) = R + A sin(i? R), where A, B are constants that 
play the role of the "amplitude" and the "frequency" of 
the oscillations respectively. However, it is far from triv- 
ial to construct and (more importantly) to theoretically 
motivate such a model, and in any case, most viable f{R) 
or other modified gravity models in the literature predict 
a smooth and very slowly evolving Gcff. 

This is the motivation for our decision to consider only 
smooth and slowly evolving models for Gcff. Therefore, 
our assumptions for the specific model for Gcff are the 
following: 

1. We assume Gcff (a) has a zero first derivative at 
a — \ and is slowly varying. 

2. We assume that Gcff (a <C 1) -^ 1, because dark 
energy should be negligible at early times. 

The second point implies that the present day Newton's 
constant has not been scaled to 1 in our conventions and 



as a result, from the moment we chose the normalization, 
any deviation of Gcff from 1 is a GR violation. Under 
these assumptions, we can approximate Gcff by a series 
expansion as follows 



Geff(a) = Geff(l) + iG;'ff(l)(a-l)2 



50 



252(a-l)' 



(3.18) 



where we have set go — GcfT(l) and 172 — Ggjj(l). Note 
that the case (30,32) — (IjO) corresponds to GR. Under 
the assumption Gcff (a <^ 1) — > 1 we can fix go = 1 ~ ^ 
and as a result we have, 

Gcff(a) = l-| + ig2(a-l)2 (3.19) 

For the choice of Eq. p. 191) we can see that the ratio 
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Whether the ratio is larger or smaller than 1 depends on 
both the sign of 32 and the sign of the integral in the last 
line of Eq. p.20|) . Since g{a) > for < a < 1, we can 
see that both the constant Iq and the integral in the last 
line of Eq. p.20p will be positive. So, when 32 < the 
ratio r of Eq. (|3.20p will be larger than 1, thus giving the 
impression that f^m qr is overestimated compared to the 
real value Jim, real and vice versa when 32 > 0. 

Finally, using the fact that gi ^ 1 we can get an ex- 
pression for 32 given a value for ilm 
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1 
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where we of course assume, based on the slowly varying 
nature of Gcff, that for the explicit example of Eq. p.l9p 
the higher derivatives gn for n > 3 are much smaller than 

52- 

Therefore, we have two options: 

• we can use Eqs. (|3.14p and (|3.17p to calculate 
^m.GR and the diagnostic respectively, 

• or we can assume some value for r2,„ and use 
Eq. (|3.2ip to calculate 32 and Gcff(<o): under the 
assumption of a slowly varying Gcff (a). 

We will explore both of these two possibilities in the next 
section. 



IV. RESULTS 

A. Tests and simulations 

In order to examine the accuracy of our method we 
performed two tests, one with simulated noise-free data 



TABLE I: Results from the tests of our analysis method for simulated noise-free data for three different values of 32 = (0, —1,1). 
In all cases the input model was Qm ~ 0.27 and as = 0.8, while Case 1 corresponds to (wojWa) = (—0.95,0) and Case 2 to 
{wo, Wa) = (—0.95, 0.5). These parameters were chosen in order to demonstrate the applicability of the method even in scenarios 
with an evolving dark energy equation of state w{a). We then fit this data to four different models for the behavior of /s8(a), 
and the numbers in brackets under each model show the recovered values for (nm,GR, 52) while the first column labeled "Original 
52" shows the assumed value of g2 by which the data were constructed. 



Original (?2 Model 1 



Model 2 



Model 3 



Case 1 



Case 2 



(0.27,0.00) 
(0.37,-0.86) 
(0.17,0.92) 



(0.27,-0.01) 
(0.37,-0.86) 
(0.17,0.91) 



(0.27,0.00) 
(0.40,-1.08) 
(0.16,1.02) 



(0.28,-0.09) 
(0.38,-0.98) 
(0.17,0.88) 



(0.28,-0.10) 
(0.38,-0.98) 
(0.17,0.87) 



(0.27,0.00) 
(0.39,-1.04) 
(0.16,0.99) 



Model 4 



(0.24,0.27) 
(0.37,-0.84) 
(0.11,1.39) 



(0.25,0.22) 
(0.38,-0.91) 
(0.12,1.35) 



and one with simulated noisy data. 

To generate the simulated noise-free data we numeri- 
cally solved Eq. (|l.ip for six models, which we divide into 
two cases depending on whether they have a constant or 
varying equation of state of dark energy and whether 
they have a constant or varying gravitational parameter: 

1. Case 1: Constant equation of state w with param- 
eters {wo, Wa) = (—0.95, 0) and three values for the 
parameter (72 ~ (—1,0, 1) 

2. Case 2: Evolving equation of state {wo,Wa) = 
(—0.95,0.5) and three values for the parameter 
52 = (-1,0,1) 

where the parameters (wo, Wa) correspond to a model for 
the dark energy equation of state w{a) = wo + Wa{l~a). 
We chose to use the values {wo,Wa) — (—0.95,0) and 
(wojWa) — (—0.95, 0.5) in our simulations in order to test 
whether our method would be sensitive to models with 
a varying equation of state of dark energy, and to ensure 
that a varying equation of state does not confuse any de- 
tection of a variation in GcS- The evolution of /s8(a) for 
each Case is shown in Fig. [TJ while in Fig. [5] we show the 
difference in fs8{a) between each Case and the ACDM 
model for f^m = 0.27 and as = 0.8. Figs.[TJi and Figs. [2^ 
correspond to Case 1, while Figs. [T]d and Figs.I^b corre- 
spond to Case 2. As it can be seen in Fig. [^b (black solid 
line), the effect of an evolving equation of state is more 
prominent, ie the departure from the expected ACDM 
model is larger, at intermediate values of the scale factor 
a, or equivalently at z ~ 1. Furthermore, as is clearly 
seen in Fig. [21 a value for 32 which is either negative or 
positive (dashed and dotted lines) causes fs8{a) to depart 
significantly at late times from its expected ACDM value. 
That is why we are able to use WiggleZ measurements of 
/s8(a) to detect non-zero values of 52 and thus Goff. 



Having established some simulated (noise-free) data, 
we then attempted to recover the values of ilm and Geff 
in that data by fitting four different models: 

1. Model 1 (Ml): The analytical solution of Eq. P3|) 
with free parameters (ilm, w). This model assumes 
GR is valid, and GcS — 1- 

2. Model 2 (M2): The generic approximation of 
Eq. (|1.7p with free parameters (J7i„,i(;) but 7 fixed 
to the approximation that is valid for constant 
equation of state 7 — ^ Z^ ■ 

3. Model 3 (M3): The generic approximation of 
Eq. ()1.7p with free parameters (51ni,7,w). This is 
the most flexible model and is appropriate for mod- 
els with varying w{a) and Gcff(a). 

4. Model 4 (M4): A series expansion fs8{a) = cio + 
020? ■ This will serve as a theory-independent test 
of the method. 

After fitting the four different models to the data by 
X^-squared minimization, we used Eqs. p.5p and (I3.2ip 
to determine the values of rim^oR and 172 respectively, 
which we show in Table HI The numbers in brackets un- 
der each model show the recovered values for (rim.GRi ^2) 
while the first column labeled "Original (72" shows the as- 
sumed value of gi by which the data were constructed. 

As it can be seen in Table HI M-^ is in all cases the most 
accurate in predicting the correct values for 92 even in the 
case of an evolving equation of state w(z). The magni- 
tude of systematic error of 52 ranges from 8% in Case 1 
(for gi ~ —1) and less than 4% for all other cases and 
choices of g^ ■ The reason for this is that M3 is a generic 
approximation to the growth factor and has enough flex- 
ibility to fit a variety of different behaviors, not only an 
evolving equation of state such as w(z) but also the more 




FIG. 1: The evolution of /s8(a) for the two different cases considered in the text. Fig. [T^ shows Case 1 that has {wo,Wa) ~ 
(—0.95,0) and 52 = (0,-1,1) while Fig. [ijj shows Case 2 that has {wo,Wa) = (—0.95,0.5). In both figures the values 
32 = (0, —1, 1) correspond to the black solid, dashed and dotted lines respectively. 
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FIG. 2: The difference in the evolution of /s8(a) between each case and the ACDM model for Qm = 0.27 and as = 0.8. Fig. [2^ 
shows Case 1 that has {wo, Wa) ~ (—0.95,0) and g2 = (0,-1,1) while Fig. [2}3 shows Case 2 that has (woyWa) = (—0.95,0.5). 
In both figures the values (72 = (0, —1, 1) correspond to the black solid, dashed and dotted lines respectively. 



complicated Case 2 that has an evolving Newton's con- 
stant (52 — (^l,!))- Surprisingly, Mi, which is the an- 
alytical solution of the growth factor for a universe with 
a constant dark energy equation of state w{a) = wq, and 
its approximation M2, do quite well in most cases with 
their magnitude of systematic error ranging from 14% in 
Case 1 with §2 — —^ and less than 2% in Case 2 with 
f/2 = —1. The polynomial M4 fares moderately in all 
cases with a magnitude of systematic error better than 
40% in both cases when 32 = 1 but manages to determine 
f/2 to within 9% of its true value in Case 2 (for 32 = — !)• 
Finally, it is worth mentioning that the reason the best- 
fitting values in Table U differ from the input values is 
that the models being fitted are only an approximation 
to the growth histories in each case, or in other words, 
the deviations are measuring the level of systematic er- 
ror inherent in using each model. Also, some of the error 
in the recovered values of 52 and r2ni,GR comes from the 
degeneracy between the effect of evolution of Newton's 
constant G^s and the effects of an evolving dark energy. 



The second test of our method was to create synthetic 
data with noise. More specifically, we created 420 mock 
WiggleZ data sets {zi , fs8i ,(Ji), where i = 1 , . . , 4. For each 
of the four data points in the 420 synthetic data sets we 
used the redshift Zi and error ai of the real WiggleZ data, 
but /s8i was created by assuming a model with an evolv- 
ing dark energy equation of state w{a) = wq + Wa{l — a) 
and a value for 32 (just like in Case 2 of the noise-free 
test), then solving Eq. (|l.ip numerically to find fs8th{a) 
and finally adding noise, ie fs8i = fs8th{zi) + Ni. We as- 
sumed that Ni is a random number drawn from a normal 
distribution with zero mean and a standard deviation 
equal to Oi. As before, we again used the same mod- 
els {Mi,M2,M^,M4) and considered the following cases 
that correspond to models with {wo,Wa) ~ (—0.95,0.5), 
fi,„ = 0.27, CTg = 0.8 and three values of 32 = (—1,0, 1). 
After fitting the four different models to each of the 420 
data sets implementing a x^-statistic, we used Eqs. p.l4p 
and (|3.2ip to determine the values of fim.GR and 52 for 
each data set respectively. In Table HIl we present the re- 



covered average values for the parameters (r^m.GR), (52), 
(0) and their respective la errors. The Itr errors quoted 
in Table |ll] correspond to one standard deviation of the 
parameters over the 420 mock data sets. 

All four models were able to recover the correct value 
for 52 within 1 cr of the actual value used to create the 
simulated data. The most accurate model, ie the one 
that was closer to the "original" value of 52 used to create 
the synthetic data, was AI3 which gave an exact value 
when 32 = —1 and was roughly within 10% and 20% 
when g2 — and 52 = 1 respectively. The most precise 
model, ie the one with the smaller error bars, was found 
to be M2 for g2 = and 32 = —1, but it was Mi for 
g2 — 1- Surprisingly and despite its simplicity compared 
to the other models, the polynomial M4 proves to be 
equally precise when 52 = —1, but its precision drops 
significantly in the other two cases ((72 = and 52 = 1) as 
the predicted error bars are larger than that of the other 
models by a factor of roughly 1.3 and 1.5 respectively. 
Also, the magnitude of the systematic error of 32 is in all 
cases roughly within 40% and 30% in the first and the last 
two cases respectively. Therefore, we find that both the 
simulations with synthetic noise-free and noisy data lead 
to the conclusion that the model M3 provides the most 
accurate measurements of the value of §2 while having a 
precision comparably good to the other models. At this 
point it is probably worth noting that the average values 
of r^m.GR and (72 in Tables HI and ITIl agree reasonably well 
each other for each case respectively. 

Finally, as it can be seen in both Tables U and [Til this 
method seems to be overestimating the value for fim.GR 
when 52 < and underestimates it when (72 > 0. So, 
when the value of flm measured by independent tech- 
niques disagrees with the value of Jlm.OR recovered using 
our technique, then this can be interpreted as evidence 
for a varying Newton's constant. This is in agreement 
with the expected behavior of the ratio r of Eq. (|3.20p . 
Furthermore, in all cases the sign of the average value of 
over the models was consistent with the assumed value 
of g2- The diagnostic would allow for a statistically- 
significant detection of varying Gofr when 32 = i 1 given 
higher signal-to-noise growth measurements. 



B. WiggleZ data analysis 

In this section we will present the results of the method 
we described in the previous sections, when it is applied 
to the real WiggleZ data. After fitting the four mod- 
els (Mi,M2,M3,Af4) to the data, we used Eqs. (pH]) . 
p.2ip and (|3.17p to determine the values of rim.GR, 32 
and respectively, which we show in Table Hill The er- 
rors of the values of the parameters were determined by 
a Fisher matrix analysis, see Chapter 15 of Ref. 25] for 
more details on the method. 

We had to perform a Fisher matrix analysis since we 
cannot estimate the errors on the parameters with the 
simulated datasets of the previous section. The reason 



for this is the specific way these synthetic data were con- 
structed. In fact, there are two ways one can create syn- 
thetic data: 

• One can assume a theoretical model, eg have spe- 
cific choices for w{a) and Goff(a), calculate fs8 at 
the WiggleZ redshifts and add gaussian noise. We 
then use these mock data to test whether the pro- 
posed method works as expected, ie it recovers the 
original and "correct" values. This can only be used 
to determine if the method works and not to derive 
the errors of parameters for the real data, since we 
did not use the latter at all. This is what we did in 
the previous section. Since our mock data resemble 
the real ones, the errors will be as well similar. 

• The second way is to draw N data points from 
the real data with replacement, ie do a boot- 
strap Monte Carlo (see Chapters 15.6.1-15.6.2 of 
Ref. 1251). This way we get datasets in which a 
random fraction of the original points are replaced 
by duplicated original points. Then one subjects 
the synthesized data to the same method and de- 
rives the desired parameters for each mock set. The 
standard deviation of the parameters corresponds 
to the error of the parameters deduced from the 
real data. 

Clearly, these two methods test different things. In our 
case, we cannot use the second method (the Bootstrap 
Monte Carlo) to determine the errors of the desired pa- 
rameters as our real dataset is too small, so we had to 
resort to Fisher matrix analysis. 

In order to confirm the Fisher matrix analysis we con- 
sidered the jack-knife approach to error estimation, see 
Ref. 26] for an introduction, and the variance between 
several simulations, see Chapter 15.6 of Ref. ^]. We 
give more details on how we implemented these two meth- 
ods in Appendix [C] We found that the method of calcu- 
lating the variance between several simulations of mock 
data is in good agreement with the Fisher estimates of 
Table IIIII On the contrary, the estimated error-bars 
from the jack-knife method are significantly larger than 
the ones obtained by the other two methods and this is 
mainly for two reasons: 



• Our estimators Eqs. (PH|) . ((XT5|) and (|X^ are 
not smooth statistics. For example, it is well known 
that the jack-knife works very well for statistics like 
the mean but fails for the median (see Ref. J2Q] ) . 

• The size of the data-set is quite small, so this affects 
the fitting procedure (deleting one point out of four 
data-points results in a data-set which is barely big 
enough to be fitted by Model M3 that has 3 free 
parameters. 

The evolution of /s8(a) for the four best-fit models and 
their la error region (gray shaded area) along with a 
ACDM model for fim = 0.27 (green dashed line) can be 



TABLE II: The results of the simulations with noisy data. All of the cases correspond to simulated data with {wo,Wa) ~ 
(—0.95,0.5), Sim = 0.27, and as ~ 0.8. The columns under each model show the recovered average values for the parameters 
(flm,GR), (52), (0) and their respective Icr errors, while the first column labeled "Original g2" shows the assumed value of 32 
by which the synthetic data were constructed. 



Original 52 Parameters 



Ml 



M2 



Ms 



Mi 



(fim,GR> 0.382 ±0.097 0.383 ±0.097 0.388 ±0.104 0.342 ±0.105 

g-2 = -1 {g-i) -0.949 ± 0.823 -0.960 ± 0.819 -1.000 ± 0.879 -0.609 ± 0.902 

(0) -0.414 ±0.360 -0.419 ±0.359 -0.437 ± 0.385 -0.266 ± 0.391 



fl2=0 



(f^m.GR) 
(ff2> 
(0) 



0.294 ±0.074 0.296 ±0.073 0.292 ± 0.076 0.236 ±0.097 
-0.206 ±0.640 -0.218 ±0.636 -0.183 ±0.663 -0.304 ± 0.845 
-0.090 ±0.273 -0.095 ±0.271 -0.080 ± 0.283 0.127 ±0.360 



(nm,GR> 0.180 ±0.050 0.181 ±0.050 0.173 ±0.053 0.122 ±0.074 

52 = ±1 {9-2) 0.82 ±0.459 0.805 ± 0.482 0.884 ± 0.489 1.316 ±0.717 

(0) 0.334 ±0.185 0.331 ±0.187 0.360±0.197 0.546 ± 0.275 



TABLE III: The results of the analysis of the WiggleZ data for g2, r^m.GR and the diagnostic 0. When 7^ 0, then this 
indicates an evolving Gcff. The evolution of /s8(a) for the best-fit models (Mi, M2, M3, M4) and their la error region (gray 
shaded area) along with a ACDM model for Sim = 0.27 (green dashed line) can be seen in Fig. [3l In each case the first row 
is the estimated result and the error using the Fisher analysis, while the second and the third rows correspond to the error 
estimated by simulated data-sets and the jack-knife respectively (see Appendix [C] for more details). The discrepancy of the 
jack-knife errors and the ones estimated by the other two approaches is explained in the text. 
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seen in Fig. [3] The most precise measurement of 52, ie 
the one that has the smaller error bars, is that of the 
model M2 that predicts 32 — ~l-128lo siO' which offers 
no convincing evidence for non-zero .92- As expected, Mi 
predicts similar results to M2, but the value of 32 for 
M3 has slightly larger errors, making its prediction even 



more compatible with GR. 

When we fix the background cosmology to ACDM , ie 
set u; = — 1, but allow the parameter 7 to be free (case 
.^3,uj=-i in Table Hill in order to capture the effects of 
an evolving Gcs , then we get a prediction which again of- 
fers no convincing evidence for a non-zero 52- The mod- 



els Ml, M2, M3 and M3_u,=_i give on average a value 
32 — —1.14 ± 0.91 which corresponds to a value for the 
second derivative of Geff(io) = -l.l9±0.95-10-'^°h^yT-^. 
Also, their prediction for the matter density (rim^GR = 
0.40 ± 0.10) is in agreement with the one measured 
from the WMAP probe (f7,n = 0.27 ± 0.03) 0]. This 
is also in agreement with the results of the diagnostic 
= —0.486 ± 0.369, which is compatible with zero. We 
should note that while the estimate of (72 depends on 
the assumption of a slowly varying GcS on the form of 
Eq. (|3.19p . the estimate of is completely model inde- 
pendent of any assumption on either Gcff or the dark 
energy equation of state w{a). 

The theory-independent model M4, fs8{a) ~ cia + 
cia} , gives a value gi = — 0.939to;96i which is in agree- 
ment with the other models, albeit with a slightly larger 
upper error, which makes it consistent with GR. Also, 
M4 can be used to get an estimate for erg. Remember- 
ing the definition /s8(a) = f{a)as{a) — -^ a S'{a) and 

the fact that for a ^ 1 we have S'{a) ~ 1 -|- 0(a), it is 



easy to see that ci = -^, so by using Eq 



we get 



(5(1) = 1 + -^ and cts = ci -I- C2/2. Using the best-fit 
values for ci and C2 we get the result (Tg = 0.75 ± 0.08, 
which is in agreement with the corresponding WMAP 7 
value as = 0.80 ±0.03 ;3]. 

We also repeated our calculations of 52 for data that 
were marginalized over an Alcock-Paczynski scale distor- 
tion factor ^] and found that the results degraded only 
by 10%, so our conclusions remained unchanged. 

Finally, the error estimates for 32 might be somewhat 
large but the big advantage of this method is the fact 
that we have not used any other data besides the four 
data points from the WiggleZ survey. Also, it is worth 
emphasizing the point that the results we get are inde- 
pendent of the model used, in that the systematic differ- 
ences in 51m. GRj ff2 and between the models are much 
smaller than the statistical errors in these measurements. 
Therefore, the only thing that limits the accuracy of this 
method is the precision of the data and the number of 
redshift bins. We will explore this issue in what follows. 



V. OTHER DATASETS AND FORECASTS 

In this section we present the results from a joint anal- 
ysis of the growth rate measurements from the WiggleZ 
and the SDSS-II data of Ref. [H], along with the op- 
timal redshift range where new data points should be 
attempted to be measured in future surveys in order to 
achieve the lowest error in 52. 

By constraining the large-scale redshift-space distor- 
tions of the LRG SDSS-II data, Ref. [2a| presented 
two new data points f{z — 0.25)a8(^ = 0.25) — 
0.3930 ± 0.0457 and f{z = 0.37)as{z = 0.37) = 0.4328 ± 
0.0370. After fitting the four models (Mi, M2,M3,M4) 
to the joint WiggleZ and LRG SDSS-II data, we used 
Eqs. (p?T4l) . ([3?2T|) and (|3Tfl) to determine the values 
of Jlni, 32 and respectively, which we show in Table 



IIVI The errors of the parameters were determined by a 
Fisher matrix analysis, see Chapter 15 of Ref. [23| for 
more details on the method. 

The most precise measurement of 52, ie the one that 
has the smaller error bars, is that of the model M2 that 
predicts 32 = "0-363]to'5367 which is compatible with 
(72 = at the Icr level. As expected. Mi predicts similar 
results to AI2, but the value of 32 for M3 has much larger 
errors, thus making its prediction compatible with both 
GR and g2 — —I at the Icr level. When we fix the back- 
ground cosmology to ACDM , ie set w = —1, but allow 
the parameter 7 to be free (case M3^u,=_i in Table ITV)) in 
order to capture the effects of an evolving Gofr, then we 
get a prediction for (72 which is compatible with both GR 
and (72 = —1 at the la level. The models Mi, M2, M3 
and M3^u;=-i give on average a value 52 — —0.34 ± 0.65 
which corresponds to a value for the second derivative 
of Gcff (io) — —3.6 ± 6.8 • 10~^-^/i^yr~^ and is again com- 
patible with GR. Also, in this case all models are com- 
patible with = at the la level. Finally, the polyno- 
mial M4 gives a measurement of a^ = 0.77 ± 0.07, which 
is in agreement with the corresponding WMAP 7 value 
as =0.80 ±0.03 [3]. 

However, we should note that one should be careful 
when combining two different data sets as this may result 
in systematic errors. For example, this was the case with 
the "Gold06" set, which was one of the early SN la data 
sets and was made of several different subsets [23]. For 
the case at hand, one reason to exercise caution would be 
some low level of correlation due to observing some com- 
mon areas of sky or the different redshift-space distortion 
models used by the two teams. Another systematic is- 
sue might be the galaxy bias, since the two surveys use 
two different types of tracers (star-forming galaxies and 
LRGs for the WiggleZ and SDSS respectively), but both 
teams have explored this assumption in detail [22|, [2^ . 

We also tried to determine the optimal redshift range 
where new data points should be attempted to be mea- 
sured in future surveys. We noticed that due to degenera- 
cies in fs8{a) there are several sweet-spots at specific red- 
shifts, ie there are points where /s8(a) is the same regard- 
less of the values of the parameters w or 172 . In Fig. Q 
(left) we show the difference in the evolution of /s8(a), 
for values of w in the range w € [—1.285, —0.755] in steps 
of Sw = 0.053 (from top to bottom), and the ACDM 
model. The range of w corresponds to five tr's from the 
WMAP7+BA0+SN la best fit w = -1.02 ± 0.053 of 
Ref. [3|. In Fig. |4] (right) we show the difference in the 
evolution of /s8(a) with w{z) = —1 and an evolving Gcff 
for values of (72 in the range (72 G [—1.5, 1.5] and in steps 
of (5(72 = 0.2 (from top to bottom) and the ACDM model. 
Both plots were created with ilm = 0.27 and as = 0.8, 
while the red points correspond to the real WiggleZ data. 
In the plot on the left, there are two sweet-spots, one at 
a — 0.365 or z = 1.74 (dashed vertical line) and an- 
other at a = 1.035 or z = —0.034 (dotted vertical line) 
which takes place in the future (the present time a = 1 
corresponds to the dashed-dotted line). On the plot on 
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FIG. 3: The evolution of /s8(o) for the best-fit models (Mi, M2, M3, M4) and their la error region (gray shaded area) along 
with a ACDM model for fim = 0.27 (green dashed line). 

TABLE IV: The results of the joint analysis of the WiggleZ±SDSS-II LRG data for 32, fim.GR and the diagnostic 0. When 
7^ 0, then this means an evolving Gcff. The models Mi, M2, M3 and M3,u,=_i give on average a value 32 — —0.34 ± 0.65. 
The results for are consistent with GR at the la level. In each case the first row is the estimated result and the error using 
the Fisher analysis, while the second and the third rows correspond to the error estimated by simulated data-sets and the 
jack-knife respectively (see Appendix [C] for more details). The discrepancy of the jack-knife errors and the ones estimated by 
the other two approaches is explained in the text. 
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the right there is only one sweet-spot at a = 0.375 or 



1.67 (dashed vertical line). These sweet spots are 
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FIG. 4: Left: The difference in the evolution of /s8(a), for values of w in the range w G [—1.285, —0.755] in steps of 5w = 0.053 
(from top to bottom), and the ACDM model. The range of w corresponds to five cr's from the WMAP7-I-BAO-I-SN la best 
fit 10 = —1.02 ± 0.053 of Ref. [3|. Right; The difference in the evolution of fs8{a) with w{z) = —1 and an evolving GcS for 
values of g2 in the range 512 £ [—1-5, 1.5] and in steps of Sg2 = 0.2 (from top to bottom) and the ACDM model. Both plots 
were created with fim = 0.27 and cts = 0.8, while the red points correspond to the real WiggleZ data. In the plot on the left, 
there are two sweet-spots, one at a = 0.365 ox z = 1.74 (dashed vertical line) and another at a = 1.035 ox z = —0.034 (dotted 
vertical line) which takes place in the future (the present time a = \ corresponds to the dashed-dotted line). On the plot on 
the right there is only one sweet-spot at a = 0.375 or z = 1.67 (dashed vertical line). Therefore, having more points at low 
redshifts (z < 0.8) will help in the detection of an evolving GcS and its discrimination from w ^ —1. In these plots the redshift 
z = 0.8 corresponds to a ~ 0.56, which is near the data point that is furthest to the left. 



really sour spots, as observations in these regions cannot 
distinguish between different values of w and 52- It is 
therefore clear that having more points at low redshifts 
(0 < z < 0.8) will help in the detection of an evolving 
Goff and its discrimination from w j^ —1. In these plots 
the redshift z — 0.8 corresponds to a ~ 0.56, which is 
near the data point that is furthest to the left. 

In order to test whether having more points at low 
redshifts can help in the detection of an evolving Geff 
and to estimate the necessary number of new points nec- 
essary to achieve this, we created synthetic data with 
noise following the same procedure as in Section IV. We 
placed 4 of the points fixed at the WiggleZ redshifts 
z = (0.22,0.41,0.60,0.78) and we also added 1, 2, 4, 8, 
12 or 14 points randomly placed in three redshift ranges 

• Range 1: < z < 0.8 

• Range 2: 1 < z < 1.4 

• Range 3: 1 < z < 1.8 

In order to keep the analysis as simple as possible, each of 
the new points was assigned an error similar to that of the 
real data at the corresponding redshifts, ie u = 0.04 for 
z > 0.3 and a = 0.07 around z < 0.3. In total we created 
18 groups of 96 synthetic data sets each, all with initial 
conditions 32 = —1, ^m = 0.27 and the rest of the param- 
eters the same as in Section IV. We then used Eq. (I3.2ip 
to determine the values of 52 and calculate the number of 
cr's to which a non-zero (72 is detected. In Fig. ^ (left) 
we show the number of extra data points in the range 



< z < 0.8 vs the corresponding number of cr's in a de- 
tection of 52 = —1- For a 2cr detection of 52 = —1, 5 new 
points are needed in the range < z < 0.8, while for a 3ct 
detection 13 new points are needed. Future observations 
in the redshift range < z < 0.8, besides the WiggleZ 
and the LRG SDSS data, also include the Baryon Os- 
cillation Spectroscopic Survey (BOSS). Specifically, they 
study the correlation functions of massive galaxies in the 
redshift range < z < 0.7 and aim to make a precise 
determination of the scale of baryon oscillations and to 
study the growth of structure and the evolution of mas- 
sive galaxies |30| . 



In Fig. [S] (right) we show the number of extra data 
points in two ranges 1 < z < 1.4 (solid line) and 1 < 
z < 1.8 (dashed line) vs the corresponding number of 
ct's needed for a detection of 5(2 = —1. The data for 
both plots were created with the same parameters as the 
synthetic data of the first row of Table |lT] (the data with 
92 = — !)■ Having more points in the redshift range 1 < 
z < 1.8 may give at maximum an 1.6cr detection of .92 = 
— 1, while having more points in the redshift range 1 < 
z < 1.4 may give a detection up to 1.7cr's. The reason for 
this is the existence of the sour spot at z ~ 1.7 described 
in Fig. m It is therefore obvious that we need both more 
and better data points in order to possibly detect an 
evolving Gcif with any statistical significance and that 
low-redshift (z < 0.8) data is most valuable. 
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FIG. 5: Left; The number of extra data points in the range < z < 0.8 vs the corresponding number of cr's in a detection of 
(/2 = —1. For a 2(7 detection of 32 = —1, 5 new points are needed in the range < 2; < 0.8, while for a 3a detection 13 new 
points are needed. Right: The number of extra data points in two ranges 1 < z < 1.4 (sohd hne) and 1 < z < 1.8 (dashed line) 
vs the corresponding number of cr's in a detection of (72 = —1. The data for both plots were created with the same parameters 
as the synthetic data of the first row of Table HTl (the data with 32 = — !)■ Having more points in the redshift range 1 < z < 1.8 
may give at most a l.Gcr detection of (72 = —1, while having more points in the redshift range 1 < z < 1.4 may give a detection 
up to l.Tcr's. The reason for this is the existence of the sour spot at 2: ~ 1.7 described in Fig. |4] 



VI. CONCLUSIONS 

We presented a method that can be used to deter- 
mine the value of the second derivative of Newton's con- 
stant G'off(to): under the assumption of a slowly vary- 
ing G'cff(a), by using the growth rate data of the Wig- 
gleZ survey. The novelty of our approach lies in the fact 
that contrary to other methods, this one does not require 
knowledge of the expansion history of the Universe H{z), 
usually found separately by using the Type la Supernovae 
data but instead we can use the generic approximation 
of Eq. (jl.7p which, as we have shown in the previous sec- 
tion, is able to fit a great variety of different scenarios 
(evolving equation of state and/or an evolving Newton's 
constant). 

In order to demonstrate the ability of our method to 
measure the correct values of 52 we performed simula- 
tions with synthetic noise-free and noisy data. We used 
four "fiducial" models, the analytic solution to Eq. (jl.ip 
Ml , the generic approximation of Eq. (|1.7p with the pa- 
rameter 7 fixed M2 or free to vary A/3 and a polyno- 
mial expansion M4. We tested the models with data cre- 
ated for different parameters {wo,Wa) of a dark energy 
equation of state w(a) = wq + Wa{i — a), three differ- 
ent values of the second derivative of Newton's constant 

52 = Geff(io)/(Geff(io)i?o') - (-1,0,1). 

All models proved to be able to measure the correct 
value of 32 within la of its true value, despite the fact 
that we tested them against data created from an evolv- 
ing equation of state w{z). The most accurate model, 
ie the one that was closer to the "true" value of 32 used 
to create the synthetic data, was M3 which gave an ex- 
act value when (72 = —1 and was roughly within 10% 



and 20% in the 32 = and (?2 = 1 respectively. The 
most precise model, ie the one with the smaller error 
bars, was found to be M2 for 52 = and 172 = —1, but 
Ml for (72 = 1. So, we found that both the simulations 
with synthetic noise-free and noisy data lead to the con- 
clusion that the model M3 provides the most accurate 
measurements of the value of 172 while having a precision 
comparably good to the other models. 

When we applied our method to the real WiggleZ data, 
we found that the models Mi, M2, M3 and M3^u,=_i give 
on average a value 172 — —1.14 ± 0.91 which corresponds 
to a value for the second derivative of Goff(to) = — 1.19± 
0.95 • 10^^"ft-^yr^^. Also, their prediction for the matter 
density (f^m.GR — 0.40 ± 0.10) is in agreement with the 
one measured from the WMAP probe (f^m = 0.27±0.03) 
,3]. Finally, we found that the polynomial M4 gives a 
measurement of erg — 0.75 ± 0.08, which is in agreement 
with the corresponding WMAP 7 value as = 0.80 ± 0.03 
-3]. 

When we performed a joint analysis of the WiggleZ and 
the SDSS-II data of Ref. [H, the models Mi, M2, M3 
and M3_^,=_i gave on average a value 32 — —0.34 ± 0.65 
which corresponds to a value for the second derivative of 
Gcs{to) = —3.6 ± 6.8 ■ lO^^^ft-^yr"^ and is again compat- 
ible with GR. Finally, the polynomial M4 gives a mea- 
surement of (78 — 0.77 ± 0.07, which is in agreement with 
the corresponding WMAP 7 value ag = 0.80 ± 0.03 Q. 

We also presented a diagnostic 



= 



Q u Y\-\ u L 



m,GR 



fl„ 



1 



n 



m,GR 



1 



1 



3/o^n 



where fim is the value of the matter density as measured 
from other independent observations. As is easily seen 
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the value of is zero only in GR and non-zero in gen- 
eral in the presence of an evolving Goff . Furthermore, we 
foimd that the sign of is correlated to the sign of the 
second derivative Goff, ie they are both either positive or 
negative. Therefore, if we assume that the value of fim is 
independently and accurately determined by other obser- 
vations, then any statistically significant deviation of the 
quantity from zero, clearly and uniquely identifies the 
presence of an evolving GeS and consequently modified 
gravity theories. We should note that while the estimate 
of g2 depends on the assumption of a slowly varying Gcff 
on the form of Eq. p.l9p . the estimate of is model in- 
dependent of any assumption on either GeS or the dark 
energy equation of state w{a). Applying it to the real 
WiggleZ data we found = -0.486 ± 0.369, which is 
consistent with zero. 

Finally, we also determined the optimal redshift range 
where new data points should be attempted to be mea- 
sured in future surveys. We noticed that due to degen- 
eracies in /s8(a) there are sour-spots at high redshifts, ie 
there points where fs8{a) is the same regardless of the 
values of the parameters w or g2- So, in terms of adding 
new data points, we found that these should be focused 
in the range < z < 0.8. 

In order to test whether having more points at low 
redshifts can help in the detection of an evolving Gcff 
and to estimate the necessary number of new points nec- 
essary to achieve this, we created synthetic data with 
noise following the same procedure as in Section IV. 
We placed 4 of the points fixed at the WiggleZ red- 
shifts z = (0.22,0.41,0.60,0.78) and we also added 1, 
2, 4, 8, 12 or 14 points randomly placed in three red- 
shift ranges < z < 0.8, 1 < z < 1.4 and 1 < z < 1.8. 
Having more data points in any of the tested redshift 
ranges of z > 1 allowed a maximum 1.7a detection of 
g2 ~ —1, while more sticking to the lower redshift range 
of < z < 0.8 allowed a much more significant detection, 
at 3cr, of 32 = — 1- 

In conclusion, we presented a method by which growth 
of structure could detect a variation in Goff in a model- 
independent manner, but found that the current growth 
rate data is consistent with a non-evolving value of Gcff. 
It is therefore imperative to have both more and better 
data points in order to possibly detect an evolving Gcff 
with any statistical significance. This makes it tantaliz- 
ing to pursue this analysis further as both the quality 
and the number of the growth rate data will increase in 
the near future. 
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Appendix A: Derivation of Eq. (|3.ip 

The growth factor satisfies the following differential 
equation: 



S"{a) 



-^'(")^ r// \ 3 17mGoff(a) 

^ («) - o75T773¥7T72^(«) 



H{a) 



2 a^H{aflHl 





(Al) 



where H{a) = | is the Hubble parameter and Eq. (JA1[) 
has the initial conditions 5{fi) = and (5'(0) = 1 for the 
growing mode. 

Eq. (|A1|) can be rewritten in the form ^^ , [13] : 



may. 

Hi 



QH{af _ 3f^niGcff(a)(5(a) 2H{aY 5"{a) 



a Hi 



Solving Eq. (|X2]) for H{aY we get 



Hi S'{a) 



(A2) 



H{af 



a''d'{ay 



m^ 



X Gcii{x)S{x)S' {x)dx 



(A3) 

Any realistic model should satisfy \ima^Q{S' (a)) -^ 1 and 
\ima^QiH{a)^a^) ^ as H{a)'^ ~ Hlfl^^Q-^ for a - 0. 
Therefore, the LHS of Eq. (jMI is equal to H{aya^6'{ay 
and we get Eq. (|3.1|) . 



Appendix B: The grow^th factor in a matter 
dominated Universe with Gcff (a) 

In a matter dominated Universe (n,ji = 1 and H{a)^ = 
H^a'^) Eq. (|1.1|) can be solved analytically for various 
cases for Goff(a). More specifically, we will consider the 
case 



Gcff (a) 



.90 



2 ^ ' 



(Bl) 



and we will show that only when go = ^^Qil^ the growth 
factor grows as 5(a) ~ a for a ^ 1. In this case, Eq. (|l.ip 
can be solved and the growing mode for the growth factor 
8(a) can be found to be: 

8(a) = a-^+'^e-^v^^^Fl(c2;l + Cl;v/3^a) 
1 



C2 



(-l + VS^-ci) 



ci = -v/1 + 24^0 + 1252 

where 1-^1(0; c; z) is a confluent hypergeometric function 
given as the limit iFi (a; c, z) = lim(,^oo 2F1 [a, b; c; z/b), 
see Ref. [l^l for more details. 

As a self-consistency test one can check that for a ^ 1, 
the growth factor scales as 

S(a) ~ a for .go = 1 - .92/2 (B2) 

S(a) ~ ai(-i+V^5+i2g.) f^j. ^^ ^ ^ (23) 



where (jB2|) corresponds to a model that does not depart 
from GR at a ^ 1 while (jB3[) corresponds to a model 
that has Gcff(O) = 1 + f . 



14 



Appendix C: Notes on the error estimation for 
Tables III and IV. 



In this section we present some details about the jack- 
knife and the variance between several simulations ap- 
proach. 

The jack-knife is a non-parametric estimation of the 
variance or more general measures of the error, see Ref. 
[26| . Assume that we have a data set of size n that 
consists of independent and identically distributed data 
points, 



Xi, A2, 



X„ 



(CI) 



and an estimator 6 = 6{Xi, X2, ... , X„). Then as- 
sume that ^(i) = 0{Xi, X2, ... , Xi^i, Xi+i, Xn) 
is the estimate after having deleted the ith point and 
^(.) = X]r=i ^(i)/""' *^^ average of the estimates 0(i). 
Then, the jack-knife estimate of the standard deviation 
is M 



1 



E 



(%)-^V))' 



(C2) 



One caveat of the jack-knife estimate is that it fails for 
non-smooth statistics, for example while it works very 
well for the mean, it fails for the median (see Ref. [2a|)- 
As seen in Table IIIIl this happens in our case as well, 
as our estimators Eqs. (|3.14p . (I3.15P and (13.211) are not 
smooth. 



For the second approach, we created synthetic data 
sets by adding gaussian noise on the WiggleZ data, ie 
mocki = {zi,fs8i,ai), where fs8i = fs8wiggieZ,i + Ni and 
the {zi,ai) correspond to the real data. Then we calcu- 
late the value of each estimator of Eqs. p.l4p . p.lSp 
and (|3.2ip and the error corresponds to the standard de- 
viation of each sample, see Chapter 15.6 of Ref. [23 for 
more details. This is different from the simulations we 
performed in Section IV. A, as that assumed a specific 
model for the synthetic data. 
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